PREFACE 


This  report  deals  with  the  summarization  of  measured  concomitant 
variables  available  from  Phase  1 of  the  Santa  Barbara  Convective  Band 
Seeding  Test  Program  conducted  for  the  Wavy  by  North  American  Weather 
Consultants  in  the  1967-68  to  1970-71  seasons.  We  are  indebted  to 
NAWC  for  provision  of  data  tapes  for  our  analyses. 

This  report  is  preliminary  to  analyses  under  way  in  reexamination 
of  the  effects  of  cloud  seeding  on  the  cited  program  and  is  part  of  a 
larger  program  to  develop  Improved  methodology  on  the  design  and  analysis 
of  weather  modification  experimentation.  It  \ned  that  data 

reported  here  will  be  used  In  association  with  ^ ommarlzatlons  of 
precipitation  data  given  in  separate  reports.  A listing  of  other 
technical  reports  issued  under  this  contract  is  appended  to  this 
report. 


Ralph  A.  Bradley 
Principal  Investigator 
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MODIFICATION  EXPERIMENT:  II.  CONCOMITANT  VARIABLES^ 


by 

Thomas  A.  Gleeson 

Department  of  Meteorology,  Florida  State  University,  Tallahassee,  Florida 
1.  Introduction. 

This  paper  presents  tabulated  values  of  concomitant  meteorological 
variables  for  Phase  I of  the  Santa  Barbara  Convective  Band  Seeding  Test 
Program  conducted  for  the  Navy  by  North  American  Weather  Consultants  in 
the  1967-68  to  1970-71  seasons.  They  were  obtained  from  data  supplied 
through  the  courtesy  of  NAWC. 

The  importance  of  such  variables  was  stressed  by  Estoque  and 
Fernandez-Partagas  [2J  In  a study  of  precipitation  dependence  on 
synoptic-scale  conditions  and  cloud  seeding:  "We  found  that  the 

variation  of  synoptic-scale  conditions  is,  by  far,  much  more  important 
than  multiple  cloud  seeding  in  determining  precipitation.  This  finding 
indicates  that  cloud  seeding  experiments  in  Florida  which  do  not  take 
into  account  the  effect  of  the  varying  synoptic  conditions,  can  lead 
to  misleading  conclusions." 

The  synoptic-scale  parameters  they  used  are  the  prevailing  flow, 
stability,  humidity  and  vertical  wind  shear.  Only  the  last  of  these 
appeared  to  have  little  value  in  their  study. 

In  the  Santa  Barbara  project,  the  concomitant  variables  deemed 
important  are  stability,  700-mb  wind  speed  and  direction,  and  500-mb 
temperature  (the  latter  being  regarded  as  an  approximation  to  cloud- 
top  temperature).  See  Thompson  £t  al  L5J  for  further  details. 
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We  have  examined  radiosonde  data  from  this  project  and  obtained 
values  of  the  following  variables:  wind  speed  and  direction,  mixing 

ratio  (a  measure  of  humidity) , Showalter  Index  (a  measure  of  stability) , 
horizontal  gradient  of  stability,  and  horizontal  transport  of  instability. 
These  are  average  values  for  atmospheric  layers.  A set  of  such  values 
is  available  for  each  of  106  bands  in  a 107-band  total  from  Phase  I of 
the  project.  All  values  are  presented  in  Table  1 along  with  data  on 
the  variables  mentioned  in  the  previous  paragraph. 

The  next  section  has  a detailed  description  of  Table  1 and  the 
variables  represented  therein.  In  the  final  section,  there  is  a dis- 
cussion of  some  statistics  calculated  from  data  in  Table  1. 

2.  Description  of  Variables  in  Table  1. 

(Numbers  in  parentheses  refer  to  rows  in  the  table.) 

(1)  Date  on  which  band  arrived  at  seeding  site. 

(2)  Time  when  band  passage  began  at  seeding  site.  Pacific  Standard  Time  (PST). 

(3)  Time  band  passage  ended  at  site,  (PST). 

(4)  Treatment  of  band.  NS  indicates  nonseeded;  S,  seeded.  If  band  was 
seeded,  number  of  fuses  ignited  is  shown. 

(5)  Number  of  radiosonde  observations  selected  to  represent  atmospheric 
conditions  associated  with  band. 

(6)  Station  that  released  radiosonde: 

SBA  - NAWC  office  at  Santa  Barbara  Ali;port 

VBG  - Vandenberg  AFB 

392  - Point  Arguello/Vandenberg 

NSI  - San  Nicolas  Island 

NTD  - Point  Mugu. 

(7)  Time  of  radiosonde  release,  (PST). 
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(8)  liixlng  Ratio,  (grams  of  water  per  1000  grams  of  dry  air) . This 
is  a measure  of  the  water-vapor  content  in  a vertical  column  from  990- 
millibars  (mb)  pressure  to  600  mb,  or  from  the  surface  to  14,000  ft., 
approximately.  Radiosonde  values  of  temperature  and  dew  point  were 
used  to  find  mixing  ratios  with  the  aid  of  a thermodynamic  diagram 
(Skew  T-Log  P Diagram) . Values  were  obtained  at  100-mb  intervals 
approximately,  and  averaged  to  give  final  value. 

(9)  700-mb  Speed  (knots)  of  the  700-mb  (10,000  ft.)  horizontal  wind, 

^ from  Thompson  et  ^ [5]. 

(10)  Direction  (degrees  east  of  north)  of  the  700-mb  wind. 

(11)  Kean  Speed  (knots)  of  horizontal  winds  in  the  layer  from  1000  ft. 
to  14,000  ft.  elevation.  This  is  the  magnitude  of  the  average  vector 
(or  resultant)  wind,  calculated  from  available  observations  at  1000-tt. 
intervals  up  to  10,000  ft.  and  at  2000-£t.  intervals  thereafter.  If 
data  were  missing  from  one  or  two  levels,  values  in  table  are  enclosed 
in  parentheses,  as  are  values  of  other  variables  (rows  12,  16,  17,  and 
18)  that  depend  on  these  observations. 

(12)  Direction  (degrees  east  of  north)  of  average  vector  wind. 

(13)  500-mb  Temperature  (degrees  Celsius),  from  Thompson  et  ^ [5J.  This 
is  an  approximate  measure  of  cloud-top  temperature,  at  about  19,000  ft. 

(14)  Stability  Claes . from  Elliott  and  Thompson  LlJ. 

UL  - unstable,  low  convective  instability  base 
UH  - unstable,  high  convective  instability  base 
ST  - stable. 

A convective  instability  base  is  the  lowest  level  in  the  atmosphere 
where  the  rate  of  decrease  of  temperature  with  height  (as  determined  from 
a radiosonde  observation)  exceeds  the  rate  of  temperature  decrease  with 
height  that  would  occur  in  a mass  of  saturated  air  at  the  same  temperature 
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as  that  of  surrounding  air.  If  It  were  lifted.  Classes  are  listed  in 
order  of  Increasing  atmospheric  stability. 

(15)  Showalter  Index;  (degrees  Celsius):  a measure  of  stability,  deter- 
mined from  radiosonde  observations  with  the  aid  of  a thermodynamic  diagram. 
The  temperature  that  a small  mass  of  moist  air  would  have  if  it  were 
lifted  from  the  850-mb  level  (approximately  5000  ft.)  to  500  mb,  is 
subtracted  from  the  observed  temperature  at  500  mb,  and  this  difference 

is  the  Index  value.  Large  positive  (negative)  values  of  the  Index 
indicate  atmospheric  stability  (instability)  , according  to  theory  and 
to  some  practical  experience  in  the  past.  See  Showalter  [4J  for  initial 
presentation  of  Index. 

(16)  Stability  Wind  (knots)  refers  to  the  magnitude  of  the  vector  V^: 

Vs  = (V9  - V^)  - (V^  - V3) 

where  V^,  , and  are  the  horizontal  wind  vectors  at  3000,  6000,  and 

9000  ft.,  respectively,  and  (V^  - V^)  and  (V^  - V^)  are  thermal-wind 
vectors  for  the  3000-6000-ft.  layer  and  9000-6000-ft.  layer  respectively. 
According  to  fairly  reasonable  physical  assumptions,  the  so-called  thermal 
wind  blows  along  a line  parallel  to  mean  isotherms  of  an  atmospheric  layer 
with  higher  temperatures  to  the  right  as  one  faces  downstream.  The  so- 
called  stability  wind  is  the  difference  between  two  thermal  winds  and 
blows  along  a line  parallel  to  lines  of  constant  atmospheric  stability 
in  the  3000-9000-ft.  layer,  with  higher  stabilities  to  the  right  as 
one  faces  downstream.  In  this  context,  higher  (lower)  stability  prevails 
in  a layer  where  the  decrease  of  temperature  with  height  is  relatively 
small  (large).  In  principle  the  magnitude  of  varies  directly  with 
horizontal  gradient  of  stability.  See  Saucier  [3J  for  discussion  of 
this  subject. 
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Note  that  effects  of  latent  heat  of  condensation  are  not  considered 
here,  but  are  considered  In  determination  of  the  Stability  Class  and 
Shovalter  Index  (rows  14  and  15). 

Values  of  Stability  Wind  speed  (and  direction)  were  calculated  from 
radiosonde  observations . 

(17)  Direction  (degrees  east  of  north)  of  Stability  Wind.  For  example, 

241  degrees  (see  Radiosonde  2,  SBA) . Tills  value  Indicates  the  most 
(least)  stable  air  Is  to  be  found  approximately  southeast  (northwest) 
of  Santa  Barbara. 

2 

(IB)  Instability  Transport  (knots  ) refers  to  the  magnitude  of  the 
vector  cross-product,  V,  x V , where  V,  and  V are  as  defined  above. 

This  Is  a measure  of  the  time  rate  at  which  less  stable  air  Is  being 
transported  by  horizontal  winds  to  the  vicinity  of  the  radiosonde  station 
and  seeding  site.  If  this  quantity  Is  negative  the  arrival  of  more 
stable  air  Is  suggested. 

To  a first  approximation,  the  6000-ft.  wind  was  used  to  represent 
the  transporting  wind  In  the  3000-6000-f t . layer  when  calculating  the 
values  displayed  In  the  table. 

3.  Discussion. 

Means  and  standard  deviations  of  concomitant  variables  were 
calculated  from  data  In  Table  1 and  are  shown  In  Table  2.  Values  were 
obtained  separately  for  nonseeded  and  seeded  bands  to  study  differences 
between  the  two  categories. 

In  a preliminary  way,  we  have  Introduced,  as  another  concomitant 
variable  In  Table  2 , the  time  of  day  when  each  band  arrived  at  the  seeding 
site.  The  reason  for  doing  so  Is  that  precipitation  tends  to  be  more 
abundant  near  noon  than  during  the  night  because  of  greater  Instability 
that  Is  produced  by  solar  heating  of  air  In  the  lowest  layers . Accord- 
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ingly,  the  time  when  band  passage  began  (Row  2,  Table  1)  was  subtracted 
from  12  (representing  noon).  The  absolute  magnitude  of  this  difference, 
expressed  to  the  nearest  tenth  of  an  hour , is  a variable  whose  summary 
statistics  are  shown  on  the  row  labelled  'Time  Difference'  in  Table  2. 

No  data  for  'Stability  Class'  are  given  in  this  table  because  the 
classes  UL,  UH,  and  ST  are  merely  descriptive.  However  one  would  expect 
some  relationship  between  these  classes  and  the  Showalter  Index  (which 
is  represented  in  Table  2)  because  both  variables  are  concerned  with 
stability.  In  Table  1,  l)L,  UH,  and  ST  appear  73,  33,  and  10  times, 
respectively,  llie  average  values  of  the  Showalter  Index  for  the  three 
groups  are  3.0,  3.3,  and  7.5,  respectively.  Indicating  a positive  re- 
lationship between  the  two  variables. 

Examination  able  2 reveals  no  major  differences  between  results 
for  nonseeded  'ed  bands.  For  each  and  every  variable  the  difference 

of  means  is  rably  smaller  than  either  of  the  standard  deviations. 

Altnough  these  differences  are  small,  they  do  exhibit  some  con- 
sistency: 

The  time  difference  between  noon  and  the  beginning  of  band  passage 
at  the  seeding  site  (Row  2)  has  a smaller  average  for  seeded  bands  (5.7) 
than  for  nonseeded  bands  (6.6). 

The  water-vapor  content  of  the  air  (Row  8)  has  a slightly  larger 
average  for  seeded  bands  (5.1  versus  5.0). 

The  mean  speed  of  the  winds  (lOGO  to  14,000  ft.),  bringing  moist 
air  inland  (Row  11),  is  greater  for  seeded  bands  (28  versus  26). 

The  Showalter  Index  (Row  15)  Indicates  less  stability  for  seeded 
bands  (3  versus  4). 

Instability  transport  into  the  region  of  study  (Row  18)  has  a larger 
average  value  for  seeded  bands  (116  versus  72). 


Taken  separately  the  effect  of  each  of  these  differences  may  be 
insignificant , but  in  combination  they  suggest  that  the  total  precipi- 
tation that  might  have  been  realized  from  seeded  bands , had  they  not 
been  seeded,  would  have  been  larger  than  the  total  amount  that  fell 
from  nonseeded  bands.  Further  analysis  of  relationships  between 
concomitant  variables  and  observed  precipitation  may  be  Informative  on 
this  issue. 

REFERENCES 

[IJ  Elliott,  R.  D. , and  J.  R.  Thompson  (1971):  Santa  Barbara  Convective 

Seeding  Test  Program  1970-71  season  and  1967-71  summary,  WWC 
Technical  Publication  5308. 

[2]  Estoque,  K.  A.,  and  J.  J.  Femandez-Partagas  (1974):  Precipitation 

dependence  on  synoptic-scale  conditions  and  cloud  seeding, 
Geoflsica  Internaclonal , 14 , 181-206. 

[3J  Saucier,  W.  J. , (1955):  Principles  of  Meteorological  Analysis. 

University  of  Chicago  Press,  438  pp.  See  p.  397. 

[4j  Showalter,  A.  K. , (1953):  A stability  index  for  thunderstorm  fore- 

casting, Bull.  Amer.  Meteor.  Soc.  . 34 , 250-252. 

[5]  Thompson,  J.  R. , K.  J.  brown,  and  R.  D.  Elliott  (1975)  : Santa 

Barbara  Convective  Band  Seeding  Test  Program  final  report. 


NWC  Technical  Publication  5804. 


OTHER  TECHNICAL  REPORTS  ON  THIS  CONTRACT 


Report  Numbers 
M-388,  ONR-110 


M-409,  ONR-111 


M-410,  ONR-112 


K-417,  ONR-117 


K-419,  ONR  118 


H-420,  ONR  119 


Hanson,  ilorgan,  bach,  Charles  L.  and  Cooley,  Edward  A. , 
Bibliography  of  Statistical  and  tieteorological  ilethod- 
ology  in  Weather  Modification,  September,  1^76. 

Bradley,  Ralph  A.  & Srivastava,  S.  S. , Correlation  In 
Polynomial  Regression,  March,  1977. 

Bach,  Charles  L. , An  Interpretive  History  of  30-years 
(1945-1975)  of  Weather  Modification,  March,  1977. 

Bradley,  Ralph  A.,  Srivastava,  S.  S.  and  Lanzdorf, 

Adolf,  Sumiiiarization  of  Precipitation  Data  In  a Weather 
Modification  Experiment:  I.  A Response  Surface  Approach, 
June,  1977. 

Gleeson,  T.  A. , Data  Sunmarlzation  In  a Weather  Modi- 
fication Experiment:  II.  Concomitant  Variables,  June,  1977. 
Hanson,  Morgan  A.  , Rank  Tests  in  'Weather  Modification 
Experiments,  June,  1977. 


-8- 


TABLE  1.  Chronology  of  Bands  and  Associated 


Instability  Transport  13t  b -106  16b  40  If  f;  58  151 


Band  Munber 


Note.  Values  in  parentheses  here  and  hereafter  .are  based  on  wind 
data  with  observations  .nissinp  fror.  one  or  two  levels. 


'?and  lunber 


23G  231  167  199  239  229  209  226  2C1  219 


142  -129  150  -103  358  430  2£6  C61  28£ 


Banl  'dumber 


110  72  -13  312  -3f7  267  -329  -117  -532  373 


aluoG  witi'.  cstcricks  were  obtained  from  VLG 
ounding  no.  17L,  released  1515  PSV. 
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;or  n Values  of  each  of  Twelve  Variables  froLi  Table  1, 
ibtained  Separately  for  Nonseeded  and  Gttc'ed  Bands  (Rows 
are  numbered  v/ith  values  corn-'^ponding,  to  rows  in 
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